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A hole in a 4d orbital of atomic xenon relaxes through Auger decay after a lifetime of 3 fs. Adding 
electronegative fluorine ligands to form xenon fluoride molecules, results in withdrawal of valence- 
electron density from Xe. Thus, within the one-center picture of Auger decay, a lowered Xe4d Auger 
width would be expected, in contradiction, however, with experiment. Employing extensive ab initio 
calculations within the framework of many-body Green's functions, we determine all available decay 
channels in XeF„ and characterize these channels by means of a two-hole population analysis. We 
derive a relation between two-hole population numbers and partial Auger widths. On this basis, 
interatomic electronic decay processes are demonstrated to be so strong in the xenon fluorides that 
they overcompensate the reduction in intra-atomic Auger width and lead to the experimentally 
observed trend. The nature of the relevant processes is discussed. These processes presumably 
underlie Auger decay in a variety of systems. 



I. INTRODUCTION 

The Auger effec1ii*2ii£iatti provides a magnificent means 
to study atoms, molecules, and surfaces. It is caused 
by a special type of electronic resonance, i.e., the decay 
by electron emission of core-ionized atoms, molecules, or 
solids. Since its discovery it has received much attention, 
because it is a fundamental process yielding deep insights 
into complex many-body effects in matter. Moreover, the 
Auger effect has proven to be useful in many experimen- 
tal situationsS*^*^ and can be used as an analytical 
tool. Therefore, an in-depth understanding of the Auger 
effect is important both for fundamental and practical 
reasons. 

A particularly interesting situation arises when sys- 
tems consisting of more than one atom are considered. 
Can the Auger decay rate of a core hole in an atom be 
influenced — or maybe even adjusted at will — by modify- 
ing the chemical environment of that atom? If yes, what 
are the underlying mechanisms? Early experiments on 
gas-phase molecules 7 and ionic solids 8 hinted at a simple 
correlation between the binding energy of a core electron 
and the Auger width of the associated core hole. For in- 
stance, withdrawal of valence electrons by electronegative 
ligands induces a chemical shift of the core level toward 
higher binding energies; at the same time, a smaller num- 
ber of valence electrons is available for Auger decay, if a 
local, atomic mechanism is assumed. Hence, within the 
one-center picture of Auger decajai, an increase in bind- 
ing energy is accompanied by a decrease in Auger width. 

Citrin 10 , however, analyzed x-ray photoemission mea- 
surements in a number of metal oxides and halides, 
and found convincing evidence for interatomic processes, 
which can lead to a trend in contradiction to the one- 
center model. Citrin's work stimulated theoretical work 
on interatomic Auger decay in crystalline NaClii and 
NaFiiS Matthew and Komninosi^ were the first to 
present a formal examination of interatomic Auger pro- 
cesses. Using strong approximations^ they concluded 
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Number of Fluorine Atoms n in XeF n 

FIG. 1: Experimental widths of the Xe4d lines in Xe, XeF2, 
XeF4, and XeF^-. The figure— suggests an average increase 
of the linewidth with an increasing number of fluorine atoms. 
Note that due to spin-orbit coupling and ligand-field effects, 
the Xe4d orbitals are not all equivalent. 



that these transitions have a small impact on the Auger 
rate, except in low-energy Auger processes. In fact, 
recent theoretical and experimental investigations show 
that in the low-energy case the effect of the chemical bond 
can be dramatici 14 i 15 i 16 i 17 i 18 

In this work, we argue that interatomic electronic de- 
cay processes can have an even greater impact than pre- 
viously assumed. We investigate as an example Auger 
decay of a Xe4d hole in Xe, XeF2, XeF4, and XeFg. It is 
demonstrated that the Auger width in the series XeF„, 
n = 0,2,4,6, would decrease drastically if only intra- 
atomic decay played a role. By contrast, the experimen- 
tally observed linewidthsi^ display — in clear contradic- 
tion to the one-center picture mentioned previously — a 
monotonic increase from Xe to XeFg. This is illustrated 
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in Fig.^ 20 We show that interatomic Auger processes are 
responsible for the discrepancy. Vibrational broadening, 
which could cause a similar effect, is very likely to be 
insignificant. 19 Single ionization spectra of the xenon flu- 
orides, which we have calculated recently>2ii2 2 i provided a 
first theoretical indication of an entirely electronic mech- 
anism. 

The xenon fluorides constitute, from our point of view, 
an ideal example, because of the relatively weak bond- 
ing between the central noble-gas atom and the fluo- 
rine ligands. This allows us to identify various inter- 
atomic decay processes and to distinguish them from 
the intra-atomic one. In particular, the purely inter- 
atomic/intermolecular electronic decay mechanisms that 
govern the inner-valence physics in weakly bound clus- 
ters of atoms or molecule o 23 ' 24 i 2 Ei 2 9i 2 7i 28 i 29 provide a con- 
ceptual framework for our analysis. Therefore, a short 
summary of these processes is given in Sec. [H] 

The nature of Auger decay in a molecule can be as- 
sessed most clearly by characterizing the available de- 
cay channels. To this end, Sec. II I II introduces a many- 
body Green's function method — the algebraic diagram- 
matic construction schem o 1 ^' 3 ? — for calculating double 
ionization spectra of molecules at a highly correlated 
level. A two-hole population analysis 1 ^! is discussed 
that serves as an important tool for determining the lo- 
calization properties of the two holes in an Auger decay 
channel. Using Wigner-Weisskopf theory^SiS we derive a 
relationship between final-state population numbers and 
partial Auger widths. The calculated double ionization 
spectra are analyzed in Sec. IIVI and the relevance of in- 
teratomic decay processes for the decay of Xe Ad vacan- 
cies is elucidated. Concluding remarks are the subject of 
Sec. E| 



II. ELECTRONIC DECAY PROCESSES IN 
WEAKLY BOUND SYSTEMS 

In this section, we review the knowledge obtained in 
the study of the electronic decay of inner-valence ionized 
clusters. The decay processes discovered in clusters, and 
the associated terminology we use in order to emphasize 
the nature of the underlying mechanisms, allow a suc- 
cinct classification of the Auger processes in XeF n . We 
mention that the decay processes discussed in the fol- 
lowing for singly ionized systems can be generalized to 
describe the decay of inner- valence vacancies of multiply 
ionized clusters 

Intra-atomic electronic decay (IAED) is usually the 
most prominent electronic decay process of core-ionized 
atoms embedded in molecules or clusters. IAED is, of 
course, the only Auger process that can take place in an 
isolated atom. In atoms, this mechanism is well under- 
stoodi 1 ! 2 ! 3 ! 4 ! 5 The initial core hole is filled by a (valence) 
electron and the excess energy is transferred via Coulomb 
interaction to a second (valence) electron, which is emit- 
ted subsequently [Fig. |2Ja)] . In molecules or clusters, 



intra-atomic decay is not the only decay process, because 
other processes can occur that involve neighboring atoms 
or molecules. 

If one removes an inner-valence electron from an iso- 
lated atom or molecule, the ionized system cannot, in 
general, decay by electron emission, for it lies energeti- 
cally below the double ionization threshold. The situa- 
tion changes if the monomer is part of a cluster. 

The initial inner- valence vacancy is filled by a valence 
electron of the same monomer and the excess energy 
is transferred to a second valence electron of a neigh- 
boring monomer. This electron is emitted subsequently. 
See Fig. EJb) for a schematic representation. The decay 
process is termed interatomic or intermolecular Coulom- 
bic decay (ICD) for clusters of atoms or clusters of 
molecules, respectively. It was identified theoretically in 
several weakly bound clusters: (HF)„ clusters ^S^SLSk 
the HF(H 2 0) 2 cluster^ (H 2 0) n clusters^ Ne„ clus- 
terS ) 24 i 28 i 37 and the NeAr dimer^ Further we would like 
to point out studies of electronic decay after inner- valence 
ionization of CN~— and of fluorinated carbanions and 
their acids ~ 

ICD of inner- valence holes in clusters is ultrafast, typ- 
ically ta king place on a time scale of the order of 10 fs 
(see Ref . liCJ for a review). Experimental evidence for ICD 
has been found recently in neon clusters S In the con- 
text of this paper, it is especially important to note that 
ICD is highly sensitive to cluster size: The ICD width 
in Ne n clusters, for instance, increases strongly with n, 
when a central neon atom carries the initial inner- valence 
vacancy and the surrounding atoms converge with n to 
the shape of the first coordination sphere in solid neon. 28 
Clearly, a larger number of nearest neighbors translates 
into a larger number of dicationic decay channels avail- 
able for ICD. 

The ionization spectra of neon and argon are quantita- 
tively very different from one another (Ref. |2(| and ref- 
erences therein). Neon has significantly higher lying sin- 
gle (IP) and double ionization potentials (DIP) compared 
to argon. The inner- valence Ne 2s IP is larger than the 
Ne -1 Ar -1 DIP in the NeAr dimeriSi Therefore, an initial 
Ne2s vacancy can lead to ICD. However, the Ne2s va- 
cancy in NeAr can decay in another way, too, because the 
Ne 2s IP is also larger than the NeAr -2 DIP. In this case, 
a valence electron from argon drops into the Ne 2s hole 
and the excess energy is transferred to another valence 
electron of argon. Upon electron emission, the dimer is 
left in a NeAr -2 state. This process is called electron- 
transfer mediated decay (ETMD)^& A schematic repre- 
sentation of the ETMD process is shown in Fig.^c). In 
Ref. |2(| it is shown that the contribution of ETMD to 
the total electronic decay width of an initial Ne 2s va- 
cancy in NeAr is appreciably smaller than the contribu- 
tion of ICD. In light of the vastly increasing ICD width in 
neon clusters^ we expect a similar size effect for ETMD, 
too. In the ETMD described, only two atoms are in- 
volved. It is called two-monomer ETMD (ETMD2). 

A three-monomer ETMD (ETMD3) has been sug- 
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(c) Two-monomer electron-transfer mediated decay (d) Three-monomer electron-transfer mediated decay 



FIG. 2: Principles of electronic decay processes in weakly bound systems: (a) intra-atomic electronic decay (IAED), (b) in- 
teratomic Coulombic decay (ICD), (c) two-monomer electron-transfer mediated decay (ETMD2), (d) three- monomer electron- 
transfer mediated decay (ETMD3). Common to all electronic decay processes is that the initial vacancy is filled by a valence 
electron from the same or a neighboring atom. The excess energy is transferred to a second valence electron of the same or 
a neighboring atom. Emission of this electron results in distinct dicationic charge distributions that allow one to characterize 
the underlying decay process. 



gested to occur for an initial Ne2s vacancy in NeA^i^i 
A schematic representation of this (proposed) process is 
given in Fig. Hfd) . A valence electron of one of the two 
argon atoms is transferred into the initial Ne2s hole. 
Energy is released, and a valence electron is ejected 
from the second argon atom, leaving the trimer in a 
NeAr _1 Ar _1 state. 

We would like to study a new example for the sur- 
prising impact of the ideas presented so far. To this 
end, we examine the decay of Xe4d vacancies in XeF„ 
and demonstrate how interatomic processes lead to an 
increase in Auger width. In order to pursue this issue, 
we need to introduce a few computational techniques and 
useful relations in the ensuing section. 



III. COMPUTATIONAL METHODS AND 
BASIC RELATIONS 

A. Algebraic Diagrammatic Construction and 
Two-Hole Population Analysis 

Many-body Green's functions are ideal tools for inves- 
tigating excitation spectra of molecular matter at var- 
ious stages of ionizationi 41 * 4 ^ The two-particle Green's 
function, specifically, provides direct access to double 



ionization spectra and associated eigenstates. Alge- 
braic diagrammatic construction (ADC) is an advanced 
perturbation-theoretical approximation scheme for eval- 
uating many-body propagatorsi 30 i 43 i 44 i 45 ADC(n) sums 
all Feynman diagrams up to nth order as well as cer- 
tain classes of diagrams up to infinite order. There is 
no Dyson equation for the two-particle Green's func- 
tion ^L^LiS -p Qe two-particle ADC evaluates the two- 
particle propagator directly in terms of a Hermitian 
eigenvalue problem^ The program 4 ^ we use to calculate 
double ionization potentials implements the two-particle 
ADC(2) schemed 

The two- hole population analysis, utilized in combina- 
tion with the ADC (2) program^ is a means to reveal 
the spatial localization of the two holes resulting from 
double ionizationii 4 * 3 ^ The analysis is carried out in a 
Mulliken-typeSi fashion by determining the atomic two- 
hole (2h) contributions to the dicationic eigenstates. To 
this end, a dicationic eigenstate in the ADC (2) scheme is 
first expanded in terms of 2h configurations in molecular 
orbital basis, 

\*»- 2 ) = J2(Xn)ij\ij) + ---, (1) 

ij 

where X n is the 2h part of the ADC eigenvector of the 
nth dicationic state \ ^^~ 2 }. The symbol \ ij) denotes 
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a spin-adapted 2h configuration with holes in molecular 
(spin) orbitals i and j, both of which are occupied in 
the Hartree-Fock ground state of the iV-electron system. 
Note that ij runs only over distinguishable 2h configura- 
tions. Each 2h configuration in molecular orbital basis, 
| ij }, can now be expanded^^i in terms of 2h configura- 
tions | paj ) deriving from the atomic orbital basis 

|y> = 2^,«l^>. (2) 

[IV 

Using this transformation, the total 2h weight X^X n (2h 
pole strength) of the nth dicationic eigenstate can be 
written aa 14 i 31 

X^X-n — Y^OY n — ^ ^ |^/xi/,n ^ ^ O puV^pa ^pa^n ; (3) 

pu per 

s / 

with eigenvector expansion coefficients y„ := C/X n in 
atomic orbital basis. O represents the overlap ma- 
trix between atomic 2h configurations, which satisfies 
U^OU=l. 

On choosing suitable sets of atomic basis functions A 
and B, one obtains population numbersi^i Qab,u = 
X) Q^v.n- Frequently, as is done in this work, A and B 

represent the atomic basis functions associated with two 
not necessarily distinct atoms within a molecule. In this 
sense, A and B can be identified with individual atoms. 
This provides a clear picture of the localization properties 
of the two holes. 



B. Relation between Two-Hole Population 
Numbers and Partial Decay Widths 

A useful approximation to the decay width of a reso- 
nance is given within the framework of Wigner-Weisskopf 
theory ; 32 ! 33 which is founded on time-dependent pertur- 
bation theory. Applications of Wigner-Weisskopf theory 
to the electronic decay of inner-valence-ionized clusters 
can be found in Refs. and|2(J There, single Hartree- 
Fock determinants were employed to describe both initial 
and final states. For more complicated systems like XeF„ 
and on the energy scale of core levels, however, the final 
states in particular require a more sophisticated treat- 
ment. 

Let the initial state be a discrete, singly ion- 

ized state approximating the core-hole resonance and a 
final state be a doubly ionized state | ) together 

with a decay electron of momentum k. For simplicity we 
employ the sudden approximation (Ref. l5ll and refer- 
ences therein) and express such a final state by the anti- 
symmetrized product ct | ^^~ 2 )j where ct is a fermionic 



creation operator^ The decay width is then given by 

r = 2- E E K <- 2 1 I *?~ x ) I 2 KE n% Ek) . 

n k 

(4) 

This equation contains three types of matrix elements: 
transition matrix element ( ^n~ 2 \ c^H | 'J'^ -1 ); initial- 
state energy E\ := (*i | H | ^i); and final-state energy 
E n% := (^- 2 \c i Hcl\^ l - 2 ). The operator H is the 
Hamiltonian in fixed-nuclei approximation. The 8 func- 
tion in Eq. ensures that the contribution to the to- 
tal decay width of only those accessible final states is 
summed that conserve energy in the decay. 
Using Eq. we obtain 

(**" a | c s H | <" x ) = E(^k (V I c % H | *f - 1 ) + • • • 

ij 

for the transition matrix elements in Eq. Let us 

assume that the initial state can be approximated by 
a single determinant l^f" 1 ) = c/|$^), where I indi- 
cates a core orbital and | &q ) the ./V-electron Hartree- 
Fock ground state. Let us further exploit the fact that 
excitations of core orbitals are, to a good approximation, 
negligible in all expansion terms of the dicationic final 
state | ^n~ 2 )- Tn e terms indicated in Eq. (JHJ with "■ • •" 
then vanish exactly due to the Slater-Condon rules^S 

To proceed, we express the dicationic electronic con- 
figurations in terms of atomic 2h vectors | jiv ) . The con- 
figurations | ij ) appearing in Eq. (JSJ) form an orthonor- 
malized set of functions, and by using the relation 

1*7 ) = E ( 07S )f<-»,P<r U P<rAj ) ( 6 ) 

we introduce an orthonormalized set of configura- 
tions \ijlv) referring to atomic orbitals. Notice that the 
atomic configurations in Eq. J2J) are not orthonormalized. 
Inserting Eq. © into Eq. JSJ leads to 

<^- 2 1 &sA i ) = E(^ f w (a* I k& I ) . 

JJ.V 

(7) 

where Y 7 [ := 0?UX n . This equation relates the tran- 
sition matrix element to a superposition of atomic-like 
quantities. 

We now insert the transition matrix element JJJ into 
the basic Eq. Q for the decay width, integrate over the 
momentum k of the emitted electron, and neglect inter- 
ference terms (cross terms). The result reads 

r = £r„, 

(8) 

The partial width r„ associated with the final dicationic 
state | ^n~ 2 ) is expressed as a sum over atomic 2h con- 
figurations. Each contributing term consists of a state- 
specific factor |^„„| 2 , which we have calculated earlier. 
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The other factor, 



where 



IT,,. 



2tt 



fe 



(fj,v\ 



N-l 



)?8{E,- El ), (9) 



depends on the dicationic state | ^„~ 2 ) only via the S 
function and is termed transition strength. Since we con- 
sider the decay of a core level, implying that the energy 
of the emitted electron is relatively high, this dependence 
on n is weak and can be neglected to a good approxima- 
tion. 

The state-specific quantities |^„ in | 2 are closely related 
to the occupation numbers discussed in Sec. IIII Al In 
principle, they could be used as an alternative defini- 
tion of occupation numbers in the spirit of Lowdin's 
population analysis^. In particular, if the atomic basis 
functions are well localized, each on its atom, we may 
put \Y'\ 2 = Q^ n and obtain 



r„ = ^|T^| 



(10) 



[At/ 



For a given molecule the quantities |T M „| 2 are univer- 
sal, and the only information that is state-specific is con- 
tained in the occupation numbers Q^, v ,n- We are able to 
compute the occupation numbers of the molecules stud- 
ied in this work, and hence Eq. Hl()[) is already a useful 
tool of analysis. 

We note, however, that many elements |T Mt/ | 2 can be of 
similar magnitude, e.g., if [i denotes a p-type atomic or- 
bital on atom A and v a p-type atomic orbital on atom B. 
We therefore group together in Eq. (|10f> the correspond- 
ing contributions: r„ = \Tab\ 2 Qn»,n + • • •< Since 

the \T^ V \ 2 grouped together are not identical, the average 
quantity |T/ib| 2 will depend somewhat on the dicationic 
state index n. Computing r = ^ T n amounts to averag- 

n 

ing over the slightly varying \Tab\ 2 and thus leads to a 
more stable result. 

Highly localized dicationic final-state holes, for exam- 
ple those that result from ICD in weakly bound clus- 
ters, lead to similar contributions to the overall decay 
width of the initially ionized state. By classifying the 
atomic pairs {pv) by their localization characteristics, 
we can form sets of atomic orbitals where \Tab\ 2 is only 
a slightly varying quantity. Let us denote with A the 
atom on which the initial core hole is localized. Then, 
if /j, and v denote atomic orbitals on A, IT^I 2 obviously 
corresponds to IAED (Sec.|TT|. Analogously, if fi belongs 
to atom A and v to atom B, we are dealing with ICD, 
and if neither \i nor v belong to A, T^ 2 describes either 
ETMD2 or ETMD3, depending on whether ,u and v be- 
long to the same atom B or to two atoms B and C. Thus, 
the minimal dissection of T into its various contributions, 
in the spirit of the above discussion, reads 



and 



Tiaed = 
Ticd = 

TeTMD2 = 
TeTMD3 = 

QlAED 
QlCD 

QETMD2 
<3eTMD3 



|Tiaed| QlAED , 
|?icd| 2 Qicd , 

|TetMD2| 2 QETMD2 
|TetMD3| 2 <3eTMD3 



n pEJ 

1/gB 

n p,,i/£B 
n f^es 



(12) 



(13) 



T — TiAED + TlCD + r 



ETMD2 



ETMD3 



(11) 



If necessary, further dissections of T into more detailed 
contributions are readily done. For instance, grouping 
together separately the contributions of /j, and v, both 
being atomic orbitals of s type, both of p type, or one 
of s and one of p type, etc. A distinction between the 
decay into singlet and triplet final dicationic states is 
also possible. Finally, we mention that the individual 
transition strengths, |Ticd| 2 , etc., may be similar in a 
series of related molecules. 



IV. RESULTS AND DISCUSSION 
A. Double Ionization Spectra 

The computed double ionization spectra for the xenon 
atom, the fluorine molecule, and the series XeF n , n = 
2,4,6, are shown in Fig. |3J The molecular ground- 
state geometries used in the calculations are described 
in Refs. l2ll andl22l It has to be stressed, however, that 
the double ionization spectrum of XeF^, unlike the single 
ionization spectrum in Refs. |^ and I22L has been calcu- 
lated assuming a geometry of Oh symmetry, instead of 
Cs v symmetry, in order to make the calculation feasible 
with the available computer resources. We also refer the 
reader to Refs. l2ll and[22l for details regarding the Gaus- 
sian basis sets employed. 

The ADC 2h pole strength plotted in Fig. [3] charac- 
terizes how well the dicationic states are described by 
2h configurations. The deviation of the 2h pole strength 
from unity yields the total contribution of three-hole- 
one-particle configurations to the state in question. Our 
focus in this work is on the energy regime below the 
Xe4<i single ionization threshold, which ranges from al- 
most 68 eV in atomic Xe to a little less than 76 eV in 
XeF 6 i 19 i 21 i 22 The double ionization spectra below 80 eV 
in Fig. |3 comprise ionization out of the outer and inner 
valence of Xe, F2, and XeF n . Since mean-field and corre- 
lation effects, as revealed by the ADC scheme, dominate 
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Double Ionization Potentials [eV] 

FIG. 3: (Color) Double ionization spectra of Xe, F2, XeF2, XeF4, and XeFg. The two-particle ADC(2) 2h pole strength is 
plotted on the ordinate to characterize how well the dicationic final states are described by 2h configurations. 



in the valence region, the spectra have not been corrected 
for relativistic effects, which are expected to be negligi- 
ble. 

Double ionization from the outer valence of the xenon 
atom takes place between 30 and 40 eV, and three dis- 
tinct DIPs are found in Fig. [3J The spectra of the se- 
ries XeF„, n — 2,4,6, display in this region increasingly 
dense lying DIPs, which is attributed to a growing num- 
ber of outer-valence states due to a higher derealiza- 
tion of the two final-state holes. The double ionization 
spectrum of F2 in Fig. [3] shows outer- valence ionization 
of F2 in the range between 40 and 55 eV. With an in- 
creasing number of fluorine ligands, the general trend in 
the spectra, especially at DIPs above 40 eV, is a decrease 
of the validity of the two-hole description of dicationic 
states of XeF„. As all fluorine atoms in XeF n are equiv- 
alent (the double ionization spectrum of XeF 6 is obtained 
in Oh symmetry), a dicationic state with a hole on xenon 
and a hole on a fluorine atom is in fact a linear combi- 



nation of configurations involving a hole on xenon and a 
hole on one of the n fluorine atoms in XeF„ . The coef- 
ficients of the linear combinations are determined by the 
strength of the interaction among the fluorine ligands, 
being the weakest in XeF2. 

The triple ionization threshold of the xenon atom^ 
is 66.2 eV; for the xenon fluorides, the triple ionization 
threshold is at least that low. Therefore, many of the di- 
cationic states we have calculated can decay to tricationic 
ones. Owing to the discrete basis set used, the appear- 
ance of a decaying state mimics a discretized Lorentzian 
curved As can be seen in Fig- EI such curves may indeed 
be identified for XeF„ above 50 eV. Interestingly, the 
Xe 4d lines in the single ionization spectra of xenon and 
its fluoridesi2i2ii2£ are thus above the triple ionization 
threshold. This indicates that there is some probability 
for the Xe4d holes to decay by emitting two electrons. 
According to a measurement by Kammerling el al.fil this 
probability amounts to about 20 % in atomic xenon. 
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When analyzing the DIPs, one has to keep in mind 
that essentially there are four contributions to a dica- 
tionic state that determine its character. The two holes 
in the dicationic state can be localized on one atom, ei- 
ther Xe or F, or on two atoms, Xe and F or two different 
F atoms. The one-site states have a large dicationic pop- 
ulation number on either xenon (Xe -2 ) or a single fluo- 
rine (F -2 ), and the two-site states are of either Xe _1 F -1 
or F1~ 1 F2~ 1 character. For each dicationic state ap- 
pearing in the spectra of the xenon fluorides, the data of 
the population analysis for equivalent fluorine atoms are 
summed up to yield the population numbers. In XeF2, 
for example, the Xe^Fl -1 and Xc~ 1 F2~ 1 population 
numbers are added to give a single Xe _1 F _1 contribu- 
tion for each dicationic state in the spectrum of XeF2 . 



B. One-site Population Numbers 

The one-site population numbers for the xenon fluo- 
rides are plotted in Fig. 0] For the graphical represen- 
tation in Fig. 01 and in Fig. the population numbers 
have been normalized, i.e., the sum of the contributions 
of Xe~ 2 , F~ 2 , Xe -1 F- 1 , and Fl"^" 1 character equals 
unity for each dicationic eigenstate. The spectra are com- 
pared to study the effect of the increasing number of fluo- 
rine atoms. The first states in the double ionization spec- 
tra possess DIPs of w 30 eV. The major contribution of 
the one-site population Xe -2 is at the lower energy part 
of the spectrum. Obviously, the density of states with 
a considerable F -2 population increases due to the in- 
creasing number of fluorine atoms. The overall distribu- 
tion of these states does not change much in the different 
compounds. Certain regions are visible where the states 
have a high F~ 2 population. These regions shift slightly 
to higher DIPs due to the reduced excess charge that the 
individual fluorine atoms acquire from the central xenon 
atom in the molecular ground state^i^ 

Conversely, the importance of Xe -2 contributions to 
the dicationic states is extremely reduced due to the with- 
drawal of valence electron density on the xenon atom. 
In XeF6, the contributions of Xe -2 character have nearly 
vanished. 



C. Two-site Population Numbers 

The two-site population numbers for the xenon fluo- 
rides are plotted in Fig. [SJ The impact of the increas- 
ing number of fluorine atoms is seen here as well. The 
states with a large F1 _1 F2 _1 population in the spec- 
trum of XeF2 are clearly separated into distinct groups of 
lines originating from Fl 2p~ x F2 2p- x , Fl 2p~ x F2 2s" 1 , 
and Fl 2s~ 1 F2 2s _1 populations. This can be concluded 
from a simple energy consideration. The IPs for ioniza- 
tion from molecular orbitals with F 2p character are lo- 
cated at ?s 20 eV, those for ionization from molecular 
orbitals with F 2s character are located at « 40 eViSLS 



In the spectrum of XeF2 in Fig. 03 the groups are approx- 
imately at 40, 60, and 80 eV. XeF2 is a linear molecule, 
the two fluorine atoms being separated by the central 
xenon atom. Such states with a large F1 -1 F2~ 1 popula- 
tion are termed opposite F1~ 1 F2~ 1 states. 

The situation changes in XeF4 and XeFg. There are 
adjacent and opposite F1 _1 F2 _1 states, and the clear 
separation between the groups seen in the spectrum 
of XeF2 is lost. One reason for this effect is the interac- 
tion between adjacent fluorine atoms, which is stronger 
than the interaction between opposite ones. This leads 
to a splitting of the fluorine lines, which is also observed 
in the single ionization spectraiSiSS 

The F1-F2 distance in adjacent F1~ 1 F2~ 1 states is 
considerably reduced in comparison to the F1-F2 dis- 
tance in opposite states, and the hole-hole repulsion 
energy varies accordingly. Of course, one has to com- 
pare states which have a large F1~ 1 F2~ 1 population 
that arises from the same types of orbitals. Such 
F1~ 1 F2~ 1 states are distributed over a small energy 
range in the double ionization spectrum. Corresponding 
lines in XeF2 mark the lower ends of such DIP ranges, 
due to the maximum distance between the two vacan- 
cies. Coulomb repulsion between holes on adjacent flu- 
orine ligands in XeF4 and XeFg results in shifts toward 
higher DIPs. 

States with a large Xe _1 F _1 population do not group 
like those of F1~ 1 F2~ 1 character. The density of these 
Xe _1 F _1 states also increases with an increasing num- 
ber of fluorine atoms, but they are distributed more 
uniformly over the whole spectral range. In XeF2 
and XeFg, these states are dominant with respect to pop- 
ulation number, and in XeF4 they are comparable to the 
F1- 1 F2- 1 states. The DIPs of the Xe^F" 1 states are 
not subject to a change of the hole -hole repulsion energy 
(in contrast to those of the F1~ 1 F2~ 1 states) because the 
Xe-F distance is the same for all fluorine ligands within a 
molecule (remember that the double ionization spectrum 
of XeFg was calculated assuming octahedral symmetry) . 



D. Xe4<2 Linewidth 

The electronic decay processes presented in Sec.HIIare 
characterized in Sec. IIII Bl according to their final-state 
population. Using Eqs. (|llfl . I|12|) . and . the analysis 
of the final-state population numbers allows one to de- 
termine the relative importance of IAED, ICD, ETMD2, 
and ETMD3 for the electronic decay of a Xe4d hole in 
XeF n . The partial decay width associated with a given 
decay process can be viewed, according to Eq. (|T2"|> . as a 
product of an averaged Coulomb matrix element and a 
two-hole population factor related to that process. 

We calculated the population factors, as defined by 
Eq. (|13() . by summing population numbers up to the en- 
ergy of the Xe 4d lines obtained in Ref . The results 
are collected in Table [I] We see that Qiaed is rather low 
in the xenon fluorides, especially when compared to the 
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Double Ionization Potential [eV] 

FIG. 4: (Color) One-site population numbers of the double ionization spectra of XeF2, XeF4, and XeFg. Each line shown is 
related to a dicationic state in the double ionization spectra in Fig. [3] The height of the line gives the one-site population 
number. Note that the contributions on xenon decrease drastically as the number of fluorine atoms increases. 



Compound 


QlAED 


QlCD 


QETMD2 


QETMD3 


Xe 


14.6 











XeF 2 


10.4 


43.8 


19.4 


22.7 


XeF 4 


7.1 


68.1 


34.3 


124.1 


XeF 6 


4.8 


80.2 


45.5 


288.6 



TABLE I: Two-hole population factors for intra-atomic elec- 
tronic decay (IAED) and for the interatomic decay processes 
ICD, ETMD2, and ETMD3 in xenon and its fluorides. The 
two- hole population factors are defined in Eq. I13H . 



population factors Qicd, <2etmd2, and Qetmd3- Hence, 
as far as population numbers are concerned, IAED is sup- 
pressed, and interatomic decay processes dominate the 
electronic decay of a Xe4c? hole. Also note the sensitive 
dependence of all population factors on the number of 
fluorine ligands. 

In order to arrive at a complete picture, the transition 
strengths in Eq. (|12(l have to be determined. For this pur- 
pose, we make a reasonable assumption: The transition 
strengths |Ttaed| 2 , |?icd| 2 , |?etmd2| 2 , and |T ET md3| 2 
are universal for xenon and its fluorides. This allows 
us — utilizing Eqs. illlfl and i|12[) — to take the averaged 
experimental Xe Ad linewidths, displayed in Fig. ^ and 



reproduced in Table II I II and the population factors in 
Table |T| to set up a linear system of four equations for 
the four unknown transition strengths. The matrix of 
population factors, however, is extremely ill-conditioned, 
as its L 2 condition number^ for instance, is 945, which 
is much larger than unity. This means that the solution 
of the linear system depends strongly on the accuracy of 
the population factors. 

Therefore, we first determine |Ttaed| 2 by solving 
Eqs. Ijlll) and 112(1 for the xenon atom. We then exploit 
that both |Tetmd2| 2 and |Tetmd3| 2 are based on elec- 
tron transfer between Xe and a fluorine ligand. They 
are thus expected to be of similar magnitude, so we 
set |Tetmd2| 2 = |?etmd3| 2 - The remaining linear system 
of three equations and two unknowns can be solved in a 
stable fashion by a least-squares procedure. 55 The calcu- 
lated transition strengths for Xe and XeF„ are displayed 
in Table [n] The results are in line with expectations: 
The transition strength associated with intra-atomic de- 
cay is greater by an order of magnitude than the inter- 
atomic energy-transfer strength of ICD. |Tetmd2| 2 and 
|7etmd3| 2 , which require interatomic orbital overlap for 
electron transfer, are even smaller. 

Finally, by combining the population factors from Ta- 
bleland the transition strengths from Table ITT1 the par- 




Decay Process 


|T| 2 [eV] 


IAED 


1.4 x 1(T 2 


ICD 


1.9 x 10~ 3 


ETMD2/ETMD3 


2.1 x 10" 4 



TABLE II: Transition strengths for electronic decay of a 
Xe4d hole in XeF„. 



tial widths Tiaed, Ftcd, Tetmd2, and Tetmd3 in XeF„ 
can be computed [Eq. I|12|)]. They are presented in Ta- 
ble lllll The conclusions that can be drawn from this table 
are truly remarkable. From Xe to XeFg, the contribution 
from IAED to the total Xe 4c? decay width drops mono- 
tonically, consistent with the one-center picture men- 
tioned in Sec. U However, that decrease is more than 
compensated for by interatomic decay, which is the rea- 
son why an overall increase in Xe4d width is observed 
in experiment. ICD is the most important interatomic 
decay mechanism. In XeF4, Ticd is already greater 
than Tiaed , and in XeFg , more than 50 % of the to- 
tal Xe4d width are caused by ICD. Decay via electron 
transfer is less significant. Nevertheless, while ETMD2 
is practically negligible throughout, Tetmd3 and Tiaed 
in XeFg are in fact comparable, due to the steep rise of 
the population factor Qetmd3 between XeF2 and XeFg 
(Table Q. 



V. CONCLUSION 

We have shown that interatomic electronic decay pro- 
cesses are of prime importance for understanding relax- 
ation of a Xe4d vacancy in the xenon fluorides. The 
experimentally observed trend of increasing Auger width 
in the series XeF„, n = 0,2,4,6, is an impressive con- 
sequence of the fact that interatomic Coulombic decay 
and electron-transfer mediated decay not only cancel out 
the drop in intra-atomic Auger rate: Interatomic decay 
dominates by far in the larger xenon fluorides. In XeFg, 
for example, 76 % of the total Xe 4c? lincwidth is due to 
the combined effect of ICD, ETMD2, and ETMD3. 

ICD in XeF„ takes place on an ultrashort time scale 
of less than 10 fs. This is of the same order of magnitude 
as the ICD lifetimes calculated for a 2s hole in a Ne 
atom surrounded by at least six neon monomers^ Our 
finding that, in XeFg, Ticd is only a little more than two 
times greater than Tetmd3, affords a new impulse for 
cluster and condensed-matter research. We anticipate 
ETMD3 to be a major decay mechanism for an inner- 
valence hole, for instance, in a neon atom embedded in 
an argon matrix. 

There are obvious consequences of our investigation for 
molecules bearing a similarity to XeF„. An interesting 
example is Au ger decay of a Si 2p hole in SiF4 . It was 
shown in Refs. [jja and[l6| that virtually not a single di- 
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Compound 


Tiaed [eV] 


Ticd [eV] 


Tetmd2 [eV] 


Tetmd3 [eV] 


r oxpt [eV] 


Xe 


0.20 











0.20 


XeF 2 


0.15 


0.08 


0.00 


0.00 


0.24 


XeF 4 


0.10 


0.13 


0.01 


0.03 


0.26 


XeF 6 


0.07 


0.15 


0.01 


0.06 


0.29 



TABLE III: Decomposition of the mean experimental electronic decay width, r cxp t, of the Xe4d- lines in Xe and XeF n into 
contributions from intra-atomic and interatomic decay processes. r Gxp t is taken from Fig. 



cationic (valence) eigenstatc in SiF4 is, at a significant 
level, of Si~ 2 character. The occurrence of primarily 
multi-center decay channels is termed foreign imaging. 
It is these channels to which ICD, ETMD2, and ETMD3 
provide access. In light of the present study, we may 
conclude that the calculation by LarkinS)££ based on the 
one-center model, must be expected to appreciably un- 
derestimate the Si 2p Auger width in SiF4. A recent mea- 
surement by Thomas et alw^ proves this expectation to 
be justified: The experimental Si 2p Auger width is more 
than five times greater than the prediction of the one- 
center model. The authors of Ref. [l8| also interpret this 
discrepancy in terms of interatomic electronic processes. 

Our approach to molecular Auger decay, as presented 
in this paper, lays its emphasis on the decay channels. 
Though not impossible, an actual computation of Auger 
widths in such complicated many-electron systems as the 
xenon fluorides still cannot be carried out in a very reli- 
able and accurate manner. We therefore put to use two 
highly powerful theoretical tools: the ADC scheme for 
the two-particle propagator and the two-hole population 
analysis of dicationic eigenstates. These two methods 
allow us to concentrate on the mechanisms underlying 
Auger decay by yielding detailed information about the 
nature of the available decay channels. 

The key relations and {T3J) between two- 



hole population numbers and partial Auger widths are, 
strictly speaking, only valid for molecules with pro- 
nounced charge localization characteristics. The xenon 
fluorides satisfy this criterion well. We are convinced, 
however, that the basic picture is much more general. In 
order to understand Auger decay in arbitrary molecules, 
it is usually insufficient to restrict the analysis to local 
valence electron densities. The physics may be even more 
complicated, but a qualitative decomposition into inter- 
atomic processes ICD, ETMD2, and ETMD3 is most 
likely still possible. 
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